Mechanisms for electron transport in atomic-scale 
one-dimensional wires: soliton and polaron effects 
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We study one-electron tunneling through atomic-scale one- dimensional wires in the presence of 
coherent electron-phonon (e-ph) coupling. We use a full quantum model for the e-ph interaction 
within the wire with open boundary conditions. We illustrate the mechanisms of transport in the 
context of molecular wires subject to boundary conditions imposing the presence of a soliton defect 
in the molecule. Competition between polarons and solitons in the coherent transport is examined. 
The transport mechanisms proposed are generally applicable to other one-dimensional nanoscale 
systems with strong e-ph coupling. 
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I. INTRODUCTION 



How will an electronic system conduct in the limit 
where it becomes entirely one-dimensional? This became 
an issue in the late 1980s, with the fabrication of semi- 
conductor structures that possessed only a single con- 
ducting "channel" , and which exhibited the phenomenon 
of quantized conductance M. With nanoscale fabrica- 
tion techniques, the class of effectively one-dimensional 
systems has expanded to include: atomic-scale structures 
on surfaces produced by scanning probe lithography , 
self-assembly j|, or directed assembly [|J; break junc- 
tions HD; car bon nanotubes [Q; and conjugated organic 
molecules p] . The measurement of the transport proper- 
ties of individual examples of such systems is now possi- 
ble; for example, conduction in individual molecules has 
recently been measured by scanning probes |9|,|i"o|| or bulk 
electrodes pjlj]; transport in nanotubes has been mea- 
sured using lithographed electrodes Similar mea- 
surements may also soon become possible for atomic- 
scale wires on surfaces. 

Many of these recent examples approach the ideal of 
a system whose atomic — as well as electronic — structure 
is truly one-dimensional. As such they may be expected 
to exhibit enhanced effects of correlations during elec- 
tron transport. Electron-electron correlations have con- 
sequences including the possible formation of Luttinger 
liquids |f|; on the other hand, the effects produced by 
electron-atom correlations can include the Peierls transi- 
tion | fl4[ , which opens a gap at the Fermi energy and ren- 
ders one-dimensional metallic systems semiconducting. If 
these quasi-one-dimensional systems are to have applica- 
tions in nanotechnology, it is vital to understand whether 
currents can pass through them, despite the possibility 



of a Peierls transition. In this paper we address this 
question, in the context of molecular wires; specifically, 
we consider finite-length poly-acetylene chains described 
by a fully quantum version of the Su-Schrieffer-Heeger 
(SSH) model fl§. 

We have chosen to study molecular conductors here, 
rather than any of the other types, because electron- 
lattice coupling in bulk samples is relatively well stud- 
ied Jl5| . The Peierls transition manifests itself in an al- 
ternation of single and double bonds (i.e., dimerisation) 
along the molecule's length, and bulk electron transport 
is dominated by mobile intrinsic defects formed when 
electrons or holes are injected into this dimerised struc- 
ture. For our purposes the most important defects are po- 
larons (which can be thought of as a local reduction in the 
dimerisation around an injected charge) and charged soli- 
tons (topological defects in the bond length alternation) . 
The soliton has an associated electronic state at mid-gap 
and the dimerisation changes sign, passing through zero 
in the neighbourhood of the defect |fj|. 

In order to understand molecular charge transport on 
nanometre lengthscales, we must develop a theory of co- 
herent transport that accounts for polaron and soliton 
formation. We took the first steps recently when we 
showed that the tunneling of carriers through molecu- 
lar wires at low temperature is dramatically enhanced 
by the formation of virtual polarons fll6|] . Therefore, 
previous treatments of tunneling in molecules based on 
clastic scattering significantly underestimate the conduc- 
tance \ \lt\ . However, the boundary conditions used in 
our previous work prevented solitons from forming in the 
molecule, as the dimerisation was constrained to be equal 
at both ends. In the present paper we study chains with 
an odd number of monomers; the ground state of such a 



* Present address: CEA-Saclay, DSM/DRECAM/SPCSI, Bat. 462, F-91191 Gif-sur-Yvette, France. E-mail: ness@cea.fr 
' E-mail: Andrew.Fisher@ucl.ac.uk 



1 



system possesses a single soliton in the centre |15],|T8| . We 
are thus able to study for the first time the role of the 
moving soliton, with its associated mid-gap state, and 
the competition between polarons and solitons in coher- 
ent transport. 



II. PHYSICAL MODEL 

The model for the molecular wires includes delocal- 
ized 7r-electrons interacting with quantum phonons. The 
electronic Hamiltonian is expressed in the basis of the 
one-electron eigenstates (labelled by n, m) of the refer- 
ence system. The e-ph coupling is linear in the phonon 
displacements and induces transitions between electronic 
states. The molecule Hamiltonian is 



{a x + 7Anm(a A + a\)clc„ 



X,n,m 



(i) 



where cjj (c„) creates (annihilates) an electron in the n- 
th electronic state with energy e„ and a\ (a\) creates 
(annihilates) a quantum of energy huj\ in the eigenmode 
of vibration A of the isolated molecule. The values of e„, 
to\ and 7Anm are calculated from the ground state of the 
neutral molecule (described by the SSH model Jl5| , [T6[ ) 
containing an odd number N of monomers. The vibra- 
tional eigenmodes V\ and frequencies of the molecule are 
calculated within the harmonic approximation fl9f| . 

To calculate the transport properties through the wire, 
the left and right ends of the molecule (atomic sites 
i = l and i = N) are connected to metallic leads 
via hopping integrals vl,r respectively. The leads are 
modeled as one-dimensional semi-infinite tight-binding 
chains with on-site energy €l,r and inter-site hopping 
matrix elements Pl.r (with no e-ph coupling). The 
scattering states \^f) for a single incoming charge car- 
rier are expanded inside the molecule onto the eigen- 
states \n, {n\}) = CnY\\{ a \) nx I ^/^- 01 tne non " 
interacting e-ph system . n\ is the occupation number 
of the phonon mode A and |0) is the vacuum state (rep- 
resenting the neutral ground state of the whole system, 



with a definite number of electrons in each part). The 
single added carrier can be anywhere in the system (on 
the left, the right, or in the molecule) and interacts with 
phonons only when inside the molecule. We only con- 
sider current-carrying states in which a single electron is 
added to the ground state; we expect this assumption to 
hold since the interval between electron transmission is 



much greater than the transit time |16|. 

We now briefly outline the calculation technique. 
The transport problem is solved by mapping the many 
body (one-electron/many-bosons) problem onto a single- 
electron problem with many scattering channels pO| , |T6| |. 
Each channel represents the different scattering processes 
by which the electron might exchange energy with the 
phonons. For an initial phonon distribution b = {m\} 
and an incoming electron from the left lead, the outgo- 
ing channels in the left and right leads are associated 
with energy-dependent reflection r a f,(e) and transmission 
iab(e) coefficients respectively (a = {n\} is the phonon 
distribution after scattering in each outgoing channel) . In 
the leads, the scattering states take the asympotic form 
of propagating Bloch waves with amplitudes r ao (reflec- 
tion) and t ao (transmission) . The wave vectors are given 
by the dispersion relations in each channel. For example, 
the energy e m of the incoming electron from the left lead 
is related to the wave vector k^ by the tight-binding- 
like dispersion relation ej„ = cl + 2/3l cos k^ . For the 
final energy e/j n of the electron transmitted to the right, 
one has: efi n = €r + 2/3rcos£^. One can then project 
out the leads from the problem and work in the molecu- 
lar wire subspace to solve for the value of the scattering 
state \*i>(E)). This state is obtained from propagating 
the source term \s(E)) (incoming electron from the left) 
via the effective Green's function G{E) defined in the 
molecular wire subspace: \^(E)) = G(E)\s(E)), where 
G is given by G(E) = [E - H - T, L {E) - ^r(E)Y_ \ H 
is the molecular wire Hamiltonian defined in Eq. (l|) and 
^l,r(E) are complex potentials arising from embedding 
the molecular wire spectrum into the continuum of states 
associated with the leads. Overall energy is conserved, 
so €i n and e/i„ are related by: E = 6i n + J2\ m \^ J \ — 
e/in + ^A n A^A 2 - 

The linear system = G\s) is solved for a finite size 



1 A similar expansion is used in the leads, ej, is then replaced 
by the charge creation operator d\ on site i. 

2 In this paper, we consider the limit of low temperatures: 
the initial phonon distribution is b = {m\} = {0} where all 
optic phonon modes are in the ground state (in this limit, 
the injection energy equals the total energy = E = 
efi n + X/a n ^^ VjJ ^)- It i s a ver y good approximation even at 
room temperature for all optic modes except the soliton trans- 
lation. For the soliton translation the condition kT <C Hcj 
restricts us to T « 230^. 
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basis set by truncating the phonon subspace (i.e consid- 
ering the lowest occupation numbers up to n™^ x in each 
mode). Furthermore the electron is coupled inside the 
molecule to a finite number iV p h of the most relevant 
phonon modes (see next section). A detailed analysis 
of the validity of such approximations can be found in 
Rcf. p6| . From the solution \^>) = G\s), one can obtain 
the reflection r a f, and transmission t a b coefficients for all 
the channels and hence the currents flowing through the 
wire. One can also calculate the expection values of any 
correlation functions between the electron and phonon 
degrees of freedom. 

We define the transmission probability T a b = 
\t a b(^in)\ 2 (iRsmk^(ef ln )/{(3 L smk^(£ m )) 3 . In our 
model, |t a (,| 2 is proportional to \{i = N, a\^)\ 2 = 
\(N, a\G\s}\ 2 and the source term is proportional to the 
velocity of the incoming electron |l6}| . Then, working in 
the real-space representation, the inelastic transmission 
probability can be rewritten in the following usual form 
: 

v 2 v 2 
T a b(efi n ,e in ) = 4 -^sinfcf (e m ) sin fcf (e fin ) 

PL PR 

x\(i = N\G ab (E)\* = 1)\ 2 , (2) 

where (N\G a b(E)\l) is the matrix element of the Green's 
function G taken between the left side i = 1 and the right 
side i = N of the molecule and the phonon configurations 
before (6) and after (a) scattering 4 . 

III. RESULTS 

First we describe the vibrational modes needed to ob- 
tain the relevant distortions of the molecule. We have 
shown that polaron formation inside the molecule af- 
ter charge injection is due to the coupling to the long- 
wavelength optic modes |Iq |. For odd- number chains 
where a soliton defect is present, we have to include in 
addition the modes responsible for the motion and the 
deformation of the soliton. We therefore include the soli- 
ton translation (the 'Goldstone mode'), the 'amplitude 
mode' related to the deformation of the soliton width, 
and higher-order deformation modes (for example the so- 
called third mode) Jl5|,^2| . Figure [j] shows the phonons 



considered for a molecule of length N = 99. As an illus- 
trative example, we also show in Fig. ^ the dimerisation 
dj for neutral isolated even- and odd-length molecules. 
The dimerisation is obtained, in terms of atomic displace- 
ments Uj, from the staggered difference between adjacent 
bond lengths: dj = (—iy (v,j+i—2uj+Uj-i). For even N, 
the neutral molecule is perfectly dimerized: the dimerisa- 
tion d c is constant in the middle of the chain despite the 
end effects. For odd AT, a soliton appears in the middle 
of the chain, acting as a domain wall separating the two 
domains of opposite sign of dimerisation. 

Now we turn on the lattice deformations induced by 
electron propagation in the (odd N) molecule connected 
to the electrodes 5 . The lattice deformations induced by 
the tunneling electron are given by the correlation func- 
tion between the electron density Pi — c\ci on site i 
and the displacement A\ = (a\ + a i x )y/Ti/2Muj\ of the 
mode A as = (i^ t A A P 4 )/(P i ). We also define the dis- 
tortion on site j due to the displaced modes when the 
electron is on site i as <^a'^ (•?')■ The dimerisa- 

tion pattern dj is then calculated from v$ . 

We plot on Fig. the absolute value of the dimerisa- 
tion dj from which the constant dimerisation d c and the 
end effects have been subtracted (the corresponding pat- 
tern for the isolated molecule is shown on Fig. ||). Such 
a choice permits us to represent the defects with more 
contrast. For an injection energy E = at mid-gap, the 
corresponding dimerisation is shown on Fig. |(a). The 
bright feature around the middle of the chain represents 
the soliton. In the absence of e-ph coupling, this would 
correspond to charge injection in resonance with the one- 
electron soliton level (see the transmission curve below) . 
The soliton would then be immobile, its position fixed 
in the middle of the chain and its width unchanged (i.e. 
straight vertical feature around j = 50). However due 
to the e-ph coupling, the soliton defect becomes mobile 
and its width varies while the electron progates. Similar 
behaviour is observed for other injection energies around 
mid-gap. As E approaches the valence band edge the 
mechanisms become different, as shown on Fig. |3|(b). 
In parallel with the soliton derealization and deforma- 
tion, one observes the formation of a polaron. This is 
characterized by the extra bright feature around the first 



3 T a i, contains as usual the ratio of the velocity in the out- 
going channel to that in the incoming channel. 

4 The factors v L R / [3l,r sin in Eq.(^) are related to the 
imaginary parts of the potentials E^.^. 

5 Calculations were performed for an electron coupled to 
the iVph = 5 lowest frequencies modes shown in Fig. [j] and 
i-octf = 4. Calculations were also done for other sets of param- 
eters. The same qualitative physics remain when the results 
are converged versus the basis set size. 
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diagonal in Fig. ||(b). In principle, the (virtual) po- 
laron corresponds to a local reduction of the dimerisation 
around the (tunneling) electron . However it appears 
as an increase (positive number) since \d^ \ is plotted on 
Fig. |^(b). Furthermore, the corresponding dimerisation 
is not simply the superposition of the polaron and the 
soliton — the two defects interact strongly together. For 
example, the dimerisation pattern for an electron at po- 
sition i w 25 is characteristic of a merging of both the 
polaron and soliton. Such mechanisms, leading to strong 
lattice distortions, would be expected to affect the trans- 
mission through the molecular wires — effects that we now 
consider. 

We define an effective total transmission probability 
T(E) arising from the contribution of the different out- 
going channels as T(E) = J2 a T a b(efin, e in ) 5(e fin + 
^2 x n\Tiuj\ — e in ) (see Fig. ||). Without e-ph coupling, 
the transmission curve presents the usual features: reso- 
nances at energies corresponding to the one-electron lev- 
els of the molecule with almost perfect transmission, and 
strong suppression of the transmission in the band-gap 
of the molecule for which propagation occurs by tunnel- 
ing. The soliton resonance (at mid-gap) is much narrower 
than the other resonances. This is to be expected since 
the corresponding one-electron (soliton) state is much 
more localized than the other states. With e-ph inter- 
action, we have already shown that a polaron can be 
formed in the molecule and that the soliton position and 
width are strongly modified. Such mechanisms affect the 
transmission in the following way: (i) the derealization 
and deformation of the soliton broaden but lower the res- 
onance peak at mid-gap, (ii) for larger injection energies, 
the formation of the polaron effectively reduces the ap- 
parent band-gap, and a polaron resonance peak appears 
(at E w 0.42 eV for the N = 99 wire) inside the original 
gap. This effective band-gap reduction correlatively in- 
creases the transmission in the tunneling regime {i.e. for 
injection energies in the gap). 



IV. CONCLUSION 

The results presented above demonstrate the com- 
plexity of electronic transport through one-dimensional 
atomic-scale wires. We have illustrated the mechanisms 
of electron transport in molecular wires when e-ph inter- 
actions are included. Owing to e-ph coupling, polarons 
can be formed inside the molecule. Polaron propagation 
is the main mechanism of transport through perfectly 
dimerised (semiconducting) molecules. The presence of 
a mid-gap state associated to a soliton defect in the mid- 
dle of the molecule involves different mechanisms for the 
transport. For injection energies around the mid-gap 
state, the delocalisation and deformation of the soliton 
is the main mechanism for electron transfer. With larger 



injection energies (still inside the energy band-gap of the 
molecule), a virtual polaron can be formed. The trans- 
port is then associated with more complicated mecha- 
nisms involving the interaction of both polaron and soli- 
ton. However in most cases, the effective reduction of the 
band-gap due to polaron formation, the derealization 
and deformation of the soliton and the polaron-soliton 
interaction increase the transmission in the tunneling 
regime. Although we studied a model derived from or- 
ganic molecules, there is good reason to believe that the 
e-ph coupling in other one-dimensional nanoscale wires 
is likely to give rise to similar phenomena. In the case of 
carriers injected into dangling-bond lines on the Si(001) 
surface 0, polaron states which are in many respects 
analogous to those in molecular wires are also formed 
p3| . Specifically, we expect the present results to be 
valid for any atomic-scale wire in which there is a degen- 
eracy between two different values of some order param- 
eter (the dimerisation in the molecular case, the surface 
dimer buckling in the dangling bond line) which are re- 
lated by a discrete symmetry and strongly coupled to the 
electrons. 
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FIG. 1. Optic component of the phonon modes for a 
chain length TV = 99. The optic component is defined as: 
(-l) j (Vx(j + 1) - 2Vx(j) + Vx(j - 1)). (a) is the mode asso- 
ciated with the soliton translation (Tia; a =0.020 eV), (b) with 
the deformation of the soliton width (hco b =0.114 eV), (c) is 
the so-called third mode (ft.w c =0.134 eV). (d,e,f) are other 
long-wavelength optic modes (Tiwa.e.f =0.144, 0.152, 0.158 eV). 




atomic position j 

FIG. 2. Dimerisation dj (in A) for isolated neutral 
molecules. Thin solid line: TV — 100, for even TV, dj is con- 
stant (d c ) in the middle of the chain. Solid line: TV — 99, for 
odd TV, a soliton exists in the middle of the chain separating 
two domains of dimerisation with opposite site. Dot-dashed 
line: the absolute value \dj\ from which d c and the end effects 
are substracted (TV = 99). Thin dashed line: deformation in- 
duced by adding an extra electron in the chain TV = 100; the 
local reduction of the dimerisation represents a polaron. 
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FIG. 3. Two-dimensional (atomic j/electron i positions) 
map of the dimerisation d}^ (in A) obtained from the atomic 
distortions uj for the N = 99 wire length. Here we plot the 
absolute value | substracting d c and the end effects. The 
incoming electron propagates from the left to the right, (a) 
Injection at mid-gap (E = 0.0): the bright feature represents 
the soliton whose position and width vary for the different 
electron positions i taken to calculate d$. (b) Injection at 
E — 0.41 eV close the polaron resonance peak in the trans- 
mission (Fig. ^): a virtual polaron is formed (bright feature 
across the first diagonal). The maximum change in |d'*'| is 
~ 0.03 — 0.04 A around the polaron and « 0.07 A around the 
soliton. 
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FIG. 4. Transmission probability T(E) versus the en- 
ergy E for a molecular wire of length N = 99. Dashed line: 
transmission through a rigid chain containing a soliton defect 
(see the sharp peak at E = corresponding to the soliton 
state resonance). Solid line with circles: transmission when 
the e-ph coupling is included. The delocalisation of the soli- 
ton and the polaron-soliton interaction broaden the resonance 
peak at mid-gap into a "mini-band" (shown in the inset, en- 
larged low-energy region). The polaron resonance peak ap- 
pears around E « 0.42 eV. 
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